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ABSTRACT 

The first observations to detect a population of distant galaxies directly in the sub- 
millimetre waveband have recently been made using the new Submillimetre Common- 
User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT). 
The results indicate that a large number of distant galaxies are radiating strongly 
in this waveband. Here we discuss their significance for source confusion in future 
millimetre/submillimetre-wave observations of both distant galaxies and cosmic mi- 
crowave background radiation (CMBR) anisotropies. Earlier estimates of such confu- 
sion involved significant extrapolation of the results of observations of galaxies at low 
redshifts; our new estimates do not, as they are derived from direct observations of dis- 
tant galaxies in the submillimetre waveband. The results have important consequences 
for the design and operation of existing and proposed millimetre/submillimetre-wave 
telescopes: the Planck Surveyor survey will be confusion-limited at frequencies greater 
than 350 GHz, even in the absence of Galactic dust emission; a la confusion noise 
limit of about 0.44 mJy beam" 1 is expected for the JCMT/SCUBA at a wavelength 
of 850 /im; and the subarcsecond resolution of large millimetre/submillimetre-wave 
interferometer arrays will be required in order to execute very deep galaxy surveys. 

Key words: galaxies: evolution - galaxies: formation - cosmic microwave background 
- cosmology: observations - diffuse radiation - radio continuum: galaxies 
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1 INTRODUCTION 

Smail, Ivison & Blain (1997, hereafter SIB) recently detected 
a new population of distant galaxies in a deep survey of 
gravitational lensing clusters in the submillimetre waveband 
using SCUBA on the JCMT (Holland et al. 1998). Subse- 
quent deeper integration in these fields, and observations of 
another five cluster fields confirm the results (Smail et al. 
1998; Blain et al. 1998). The identification of the bright- 
est 850- /im source with a dust-shrouded starburst/active 
galactic nucleus at redshift z — 2.8 (Ivison et al. 1998) and 
the non-correspondence of the detected sources with cluster 
galaxies support these claims. Prior to these observations the 
surface density of faint galaxies in this waveband was uncer- 
tain by up to three orders of magnitude (Blain & Longair 
1993, 1996, hereafter BL96), but it is now known to an accu- 
racy of about 50 per cent. The estimation of the magnitude 
of confusion noise (Scheuer 1957) due to the submillimetre- 
wave continuum radiation of discrete extragalactic sources 
has previously required extensive extrapolation from the 
results of observations of low-redshift galaxies (Longair & 
Sunyaev 1969; Franceschini et al. 1989; Helou & Beichman 



1990; Fischer & Lange 1993; Toffolatti et al. 1995; Gawiser 
& Smoot 1997), and so the uncertainties in the estimates ob- 
tained were large as compared with those in the radio and 
far-infrared wavebands. Our new direct observations allow 
more reliable estimates to be made. 

First, we briefly discuss our observations at 850 /im. We 
introduce a simple parametric model of galaxy evolution, 
which is similar in form to the IRAS-based models discussed 
by BL96, but fitted to the counts derived by SIB, to pre- 
dict the properties of confusing sources at other millime- 
tre/submillimetre wavelengths. These models include only 
the contribution from dusty extragalactic objects. At the 
longest and shortest wavelengths in this waveband a sig- 
nificant contribution to confusion is expected due to other 
sources of radiation - discrete extragalactic non-thermal ra- 
dio emission and Galactic dust emission respectively. Sec- 
ondly, we present our estimates of confusion noise as a func- 
tion of both the angular scale and wavelength of observa- 
tions. Finally, we discuss the relevance of confusion for obser- 
vations of both distant galaxies and the CMBR using exist- 
ing and planned millimetre/submillimetre-wave telescopes: 
the JCMT/SCUBA; the airborne Stratospheric Observatory 



© 0000 RAS 



2 A. W. Blain, R. J. Ivison & Ian Smail 

Table 1. Instrumental parameters for instruments and telescopes 
shown in Fig. 2, and a large MIA. References can be found in 
Section 1. It is proposed to attach a prototype of the BOLOCAM 
instrument destined for the LMT to the 10.4-m Caltech Submil- 
limeter Observatory (CSO). Unless otherwise stated the sensitiv- 
ities quoted refer to a la detection in a 1-h integration including 
overheads. 



Instrument 


A / 


T^A/tTl\/T Roam 
r vv nivi oeam- 


Sensitivity 




fim 


width / arcmin 




SCUBA 


450 


0.12 


50mjy 


(Jiggle-map) 


850 


0.22 


8mJy 


FIRST (25- 


250 


0.30 


3.3 mjy 


arcmin map) 


350 


0.42 


2.9 mjy 




500 


0.60 


2.4 mjy 


Planck Surveyor 


350 


4.4 


26 mjy 


(Nominal 14- 


550 


4.4 


19mJy 


month mission) 


850 


4.4 


16 mjy 






i i 
/ .1 


11 mjy 




2000 


10.3 


11 mjy 


SOFIA 


450 


1.3 


3.5 mjy 


(Point source) 








BOOMERANG 


770 


20 


AT/T = 


(Nominal 5-day 






3.6 x 10" 6 


mission) 


1250 


20 


As above 


BOLOCAM/CSO 


1100 


0.50 


0.5 mjy 


(36-arcmin 2 map) 


1400 


0.63 


0.5 mjy 




2100 


0.95 


0.5 mjy 


BOLOCAM/LMT 


1100 


0.10 


0.033 mjy 


(2-arcmin 2 map) 


1400 


0.13 


0.033 mjy 




2100 


0.20 


0.033 mjy 


MIA 


870 


0.02 


0.028 mjy 


(single pointing) 


1300 


0.02 


0.017 mjy 



For Infrared Astronomy (SOFIA: Becklin 1997); the Large 
Millimeter Telescope (LMT/GTM), a ground-based 50-m 
single-antenna telescope (Schloerb 1997) that will incorpo- 
rate the BOLOCAM bolometer array receiver (Mauskopf & 
Bock 1997); large ground-based millimetre/submillimetre- 
wave interferometer arrays (MIAs: Brown 1996; Downes 
1996); the ESA Far-Infrared and Submillimetre Telescope 
(FIRST: Pilbratt 1997) and Planck Surveyor (Bersanelli et 
al. 1996) space missions, and balloon-borne instruments, for 
example BOOMERANG (Lange et al. 1995). Parameters of 
these instruments and telescopes are listed in Table 1. An 
advanced-technology 10-m ground-based telescope is also 
planned for the South Pole (Stark et al. 1998). 



2 THE POPULATION OF GALAXIES IN THE 
SUBMILLIMETRE WAVEBAND 

The galaxies discovered by SIB were detected in the fields 
of two clusters, Abell370 and C12244— 02, each of which was 
observed for 25 ks at wavelengths of 450 and 850 /xm in a 
circular field 2.3 arcmin in diameter. The la noise flux den- 
sities in these maps were 14 and 2 mjy respectively. In the 
total field area of about 10 arcmin , six sources were de- 
tected at 850 /im and one at 450 /jm. The results presented 
in this paper are based on these published observations; how- 
ever, we have now expanded our survey to cover five addi- 
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Limiting Flux Density S / mjy 

Figure 1. Predicted counts of galaxies in the millimc- 
tre/submillimetre waveband, based on the model discussed by 
SIB. The model is in agreement with recent observations by 
Kawara et al. (1997 - K) and the la upper limit of Wilner & 
Wright (1997 - WW). Note that the errors on the SIB counts arc 
not independent. 

tional clusters, and increased the total field area to more 
than 31 arcmin 2 (Smail et al. 1998; Blain et al. 1998). The 
magnification bias due to gravitational lensing in the clus- 
ters (Blain 1997b) was exploited by SIB in order to assist 
the detection of background galaxies; however, these results 
should soon be confirmed by blank-field SCUBA surveys 
(BL96; Pearson & Rowan- Robinson 1996). After correct- 
ing for the moderate magnification bias we estimate that 
the surface density of galaxies with flux densities larger 
than 4mJy is (2.5 ± 1.4) x 10 3 deg -2 at a wavelength of 
850 /jm. This surface density is larger than that predicted 
in the most strongly evolving model by BL96, but is nev- 
ertheless consistent with both the observed intensity of dif- 
fuse submillimetre-wave background radiation (Puget et al. 
1996; Schlegel, Finkbeiner & Davis 1998; Hauser et al. 1998; 
Fixsen et al. 1998) and the abundance of heavy elements at 
the present epoch. In Fig. 1 the cumulative counts of galax- 
ies inferred from the 850-/im observations are compared with 
the well- fitting model from SIB, which is also used to predict 
the corresponding counts at the wavelengths of four other at- 
mospheric windows - 2800, 1300, 450 and 350 fim - and at a 
wavelength of 175 /jm, which is inaccessible to ground-based 
telescopes. This model is described in SIB, and assumes that 
the local population of IRA S galaxies (Saunders et al. 1990) 
undergoes pure luminosity evolution of the form (1 + zf 
out to z = 2.6, with a fixed evolution factor of 46.7 in the 
range 2.6 < z < 7. The counts predicted by this model at 
175 /im and 2.8 mm are consistent with recent observations 
by Kawara et al. (1997) and Wilner & Wright (1997): see 
Fig. 1. The predicted counts are expected to be most accu- 
rate at the observed wavelength of 850 /im. 



3 SOURCE CONFUSION 

Source confusion from faint galaxies that are not individu- 
ally recognized can be calculated rigorously using the for- 
malism derived by Scheuer (1957). However, for reasonable 
source count models confusion is expected to be most sig- 
nificant when the cumulative surface density of sources N 
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Figure 2. The 1<t confusion noise predicted using the models shown in Fig. 1 as a function of both observing frequency and angular 
scale. In (a) the noise equivalent flux density S v is plotted. In (b) the relative thermodynamic temperature uncertainty in the 
CMBR spectrum AT/T is plotted. The contours are spaced by factors of 10, but only the four lowest contours are labelled in 
(b). The resolution limits and observing frequencies of several telescopes (Table 1) are also plotted using different point styles. In 
(a) the two-digit numbers alongside the points represent the logarithms of the integration times in hours after which confusion 
noise is expected to exceed instrumental noise for each telescope; for comparison, 1 day, 30 days and 1 year correspond to values of 
1.38, 2.86 and 3.94 respectively. For BOLOCAM, BOOMERANG, FIRST, SOFIA and SCUBA these times refer to an imaging 
observation in a single pointing. The ratio of the integration times for imaging and point source photometry assumed for SOFIA 
is 20, consistent with the estimates for the FIRST bolometer camera (Griffin 1997). The time for Planck Surveyor refers to the 
noise per pixel achieved in an all-sky survey. The resolution limit of a large ground-based MIA would lie below the bottom of both 
panels. The corresponding two-digit numbers for an MIA are 0.9 at 353 GHz and 1.9 at 230 GHz, assuming a 3-arcsec beam. 



is approximately lbeam -1 ; the flux density S at this sur- 
face density is roughly equal to the la confusion noise. At 
larger surface densities many sources are expected within 
each beam, and so the point-to-point fluctuations in the de- 
tected signal due to the absence or presence of individual 
sources in the beam is reduced. At smaller surface densi- 
ties, a smaller number of beams contain a source, and so the 
significance of confusion is again reduced. 

The model responsible for the counts shown in Fig. 1 
is used to obtain estimates of confusion noise in the mil- 
limetre/submillimetre waveband due to discrete extragalac- 
tic sources of dust continuum radiation. Given that the 
observed counts are currently uncertain by 50 per cent, 
it is reasonable to identify the la source confusion noise 
with the flux density at which a count of 1 beam -1 is pre- 
dicted. A top-hat beam is assumed in the calculations, and 
so other authors can readily adapt the results to fit with 
other beamshapes. Simulations of confusion noise carried 
out by randomly placing sources drawn from the model 
counts into Gaussian beams predict values several times 
larger than the values produced by the lbeam -1 analysis. 
Hence, while this paper does not present the definitive calcu- 
lation of millimetre/submillimetre-wave confusion noise, it is 



the first to be observationally-based, makes reliable predic- 
tions that represent the best state of current knowledge and 
improves very considerably on the important work of earlier 
authors. 

The resulting estimates of confusion noise as a func- 
tion of both observing wavelength and resolution are pre- 
sented in Fig. 2. The specified angular resolutions and ob- 
serving frequencies of a range of telescopes and instruments 
- the JCMT/SCUBA, FIRST, Planck Surveyor, BOLO- 
CAM, BOOMERANG and SOFIA - are also shown and 
listed in Table 1. The integration time required for confusion 
noise to dominate the noise in an integration is also given 
for all the instruments. The power-law slopes of the counts 
in Fig. 1, from which these estimates are determined, are 
about (3 — —2, where N oc S 13 , and so confusion is expected 
to be insignificant at flux densities about five times larger 
than the confusion noise estimates presented in Fig. 2(a). 

At frequencies less than about 100 GHz on angular 
scales larger than about 0.2°, discrete non-thermal radio 
sources are expected to start to dominate the confusion 
noise. At the highest frequencies shown in Fig. 2, Galactic 
emission is likely to be the most significant source of con- 
fusion noise. The contribution of these sources to confusion 
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is not included in Fig. 2, but extrapolated estimates can be 
found in Franceschini et al. (1989) and Helou & Beichman 
(1990). 

The magnification bias owing to lensing is expected to 
modify the effects of confusion in the fields of clusters. If 
the inner 1-arcmin region of the cluster is resolved, then the 
confusion limits obtained after correcting for the magnifica- 
tion bias are expected to be reduced by a factor of about 
2 as compared with blank fields, because lensing increases 
both the flux density of background galaxies and their sepa- 
ration on the sky. However, if clusters cannot be resolved on 
arcminute scales, then the confusion limit is expected to in- 
crease by a factor of up to about 3. The increase is expected 
to be largest if the observing beam is matched to the Ein- 
stein radius of the cluster, and to be smaller in larger beams. 
The effect becomes negligible on scales about 10 times larger 
than the Einstein radius. These effects are discussed in more 
detail elsewhere (Blain 1998a). 



4 CONSEQUENCES FOR OBSERVATIONS 

4.1 Submillimetre-wave blank-field galaxy surveys 

SCUBA, which operates at wavelengths in the range 2 mm to 
350 fim, is currently the most sensitive submillimetre-wave 
instrument on any telescope. Blank-field surveys are being 
carried out using this instrument, and so it is important to 
quantify the flux density at which confusion is expected to 
contribute a significant source of noise to such observations. 
The lcr equivalent flux density of confusion noise in SCUBA 
observations at wavelengths of 450 and 850 //m is expected 
to be about 0.14 and 0.44 mjy beam -1 respectively. Hence, 
confusion will impose an effective limit to the depth of a 
blank-field 850-^im SCUBA survey, but the confusion limit 
will only be reached after about 40 8-hour shifts in a sin- 
gle 2.3-arcmin field, including overheads, and so does not 
present a problem to planned surveys. 

We predict that the la confusion limits for the bolome- 
ter instrument on the FIRST mission (Griffin 1997) arc 
about 2.4, 2.9 and 3.3 mjy at frequencies of 600, 857 and 
1200 GHz respectively. In a single 36-arcmin 2 field these 
limits should be reached in integrations lasting about 0.4, 

1.2 and 5.6 h respectively. Thus relatively shallow large-area 
galaxy surveys should be the natural goals of the FIRSTmis- 
sion. SOFIA is expected to be confused at a la flux density 
of about 5 mjy at a wavelength of 450 /im; however, SOFIA 
is expected to be less sensitive as compared with FIRST 
(Becklin 1997), and so the confusion limit would only be 
confronted in a considerably longer integration: see Fig. 2(a). 
Future ground-based interferometer arrays with subarcsec- 
ond resolution should not suffer from the effects of confusion 
by extragalactic sources unless sensitivities of a few /xJy are 
obtained. The prospects for using such arrays for extremely 
deep galaxy surveys are excellent (Blain 1996). Wide-field 
millimetre-wave surveys to a confusion limited depth of or- 
der 1 and 0.05 mjy will be possible when the BOLOCAM 
receiver is fitted to the CSO and LMT/GTM respectively. 
The large mapping speed of BOLOCAM ensures that it is 
an ideal instrument to detect sources for subsequent study 
using the fine angular resolution of a MIA. 



Table 2. The la source confusion noise (AT/T) con f expected 
from dusty galaxies in each of the five observing bands of the 
Planck Surveyor HFI instrument and the specifications of the 
instrument, including its intrinsic sensitivity (AT/T) sons . 



A / mm 


v 1 GHz 


Pixel Area 
/ arcmin 


(AT/T) sens 


(AT/T) conf 


2.00 


150 


106 


1.2xl0~ 6 


2.4xl0~ 7 


1.38 


217 


51 


2.0xl0~ 6 


l.OxlO -6 


0.85 


353 


19 


1.2xl0~ 5 


8.1 xl0~ 6 


0.55 


545 


19 


7.7xl0 -5 


8.1 xl0~ 5 


0.35 


857 


19 


4.2 x!0~ 3 


4.5xl0 -3 



4.2 Observations of the CMBR 

An all-sky millimetre/submillimetre-wave survey of fluctua- 
tions in the intensity of the CMBR is planned at an angular 
resolution of about 5 arcmin and a sensitivity AT/T ~ 10~ 6 
using the High Frequency Instrument (HFI) on the Planck 
Surveyor mission. The predicted confusion noise in each of 
the five frequency channels of the HFI is presented in Ta- 
ble 2. These estimates are comparable to the intrinsic sen- 
sitivity of the HFI in the nominal mission (Bersanelli et 
al. 1996) at the three highest frequencies (353, 545 and 
857 GHz); however, confusion noise from dusty galaxies is 
expected to make only a minor contribution to the total 
noise in the maps at the two lowest observing frequencies 
-150 and 217 GHz. The significant confusion expected at 
high frequencies will be increased by an additional factor 
due to lensing in the fields of clusters (Blain 1997a, 1998b) , 
and so Planck Surveyor may be less efficient at detecting 
the Sunyaev-Zel'dovich effect in clusters than suggested by 
earlier studies (Bersanelli et al. 1996; Haehnelt & Tegmark 
1996; Aghanim et al. 1997) 

Stratospheric balloon-borne instruments also observe 
the CMBR in the millimetre/submillimetre waveband. The 
two highest observing frequencies of the BOOMERANG 
instrument, which has a 20-arcmin beam, are 240 and 
390 GHz. la source confusion noise of AT/T = 3.0 x 10~ 7 
and 1.8 x 10 is expected at these frequencies, respectively, 
and so source confusion noise is expected to be significant 
only at 390 GHz in an extremely long-duration flight. 

There appear to be good prospects for making a ground- 
based observation of CMBR anisotropies on small angular 
scales using the lowest frequency channel of the BOLOCAM 
receiver at the CSO: see Fig. 2(b). la source confusion noise 
of AT/T — 3.9 x 10~ 6 is expected in this channel, at a 
wavelength of 2.1 mm. Observations of the same field using 
the two higher frequency channels, at wavelengths of 1.4 
and 1.1mm should allow Galactic emission and any bright 
foreground sources to be subtracted from the resulting map. 

Although millimetre/submillimetre-wave observations 
of the CMBR will be complicated by the presence of lumi- 
nous distant dusty galaxies, these sources could be of great 
interest for both the determination of the star formation rate 
in the early universe (Blain et al. 1998) and the compilation 
of large samples of gravitational lenses (Blain 1996, 1997c, 
1998a). 



5 CONCLUSIONS 
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(i) We have estimated the source confusion noise expected 
in millimetre/submillimetre-wave observations, based on di- 
rect observations of distant galaxies at a wavelength of 
850 /im (SIB). Although comparable in magnitude to pre- 
vious estimates of the confusion noise, our predictions are 
much less dependent on either extrapolations or model- 
dependent factors. 

(ii) Reliable limits to the maximum depth of galaxy sur- 
veys that can be carried out using either ground-based or 
space-borne telescopes are imposed by these estimates. For 
example, a la sensitivity limit of about 0.4 mjy is pre- 
dicted for the JCMT/SCUBA at 850 /im. This limit could 
be reached in a fully-sampled map after about 300 h of in- 
tegration in a single field, and so is unlikely to present a 
problem to planned surveys. Future airborne and space- 
borne submillimetre-wave telescopes are expected to be con- 
fusion limited much more quickly, and so extremely deep 
millimetre/submillimetre-wave galaxy surveys will require 
the subarcsecond resolving power of large ground-based in- 
terferometer arrays. 

(iii) Our estimates will be useful for planning and op- 
crating future millimetre/submillimetre-wave telescopes. A 
FIRST survey should cover a large area of sky and not in- 
volve integrations longer than a few hours in any single- 
pointed field. Confusion noise arising from distant dusty 
galaxies will dominate instrumental noise in the planned 
all-sky survey using the Planck Surveyor HFI instrument 
at frequencies greater than 350 GHz. Confusion noise will be 
more severe in the fields of clusters, reducing the efficiency of 
Planck Surveyor for the detection of the Sunyaev-Zel'dovich 
effect. 
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